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t^j- . The quantum fluctuations of horizons in Robertson- Walker universes and 

in the Schwarzschild spacetime are discussed. The source of the metric fluctu- 
ations is taken to be quantum linear perturbations of the gravitational field. 
Lightcone fluctuations arise when the retarded Green's function for a massless 
q-i field is averaged over these metric fluctuations. This averaging replaces the 

delta-function on the classical lightcone with a Gaussian function, the width 
of which is a measure of the scale of the lightcone fluctuations. Horizon fluc- 
tuations are taken to be measured in the frame of a geodesic observer falling 
through the horizon. In the case of an expanding universe, this is a comoving 
5^ ■ observer either entering or leaving the horizon of another observer. In the black 



hole case, we take this observer to be one who falls freely from rest at infin- 
ity. We find that cosmological horizon fluctuations are typically characterized 
by the Planck length. However, black hole horizon fluctuations in this model 
are much smaller than Planck dimensions for black holes whose mass exceeds 
the Planck mass. Furthermore, we find black hole horizon fluctuations which 
are sufficiently small as not to invalidate the semiclassical derivation of the 
Hawking process. 
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1 Introduction 



One of the characteristics of classical gravitation is the existence of horizons, surfaces 
which divide spacetime into causally distinct regions. The most striking example is 
the black hole horizon, the boundary which hides the events within from the outside 
world. Cosmological models also possess horizons of a different sort; a given observer 
generally cannot see all of the other observers in the universe at a given time. If the 
expansion rate in comoving time is less than linear, then previously unseen objects 
enter the observer's horizon. If it is faster than linear (inflationary expansion), then 
objects leave the horizon. Horizons are of course lightcones, and the notion of an 
event being within or without a horizon means being at a timelike or a spacelike 
separation, respectively. 

It is expected that quantum metric fluctuations should smear out this precise 
distinction, and hence smear out the classical concept of a horizon. Information 
could presumably leak across the horizon in a way that is not allowed by classical 
physics. Bekenstein and Mukhanov |IJ have suggested that horizon fluctuations could 
lead to discreteness of the spectrum of black holes. Several other authors 0, [| have 
recently made proposals for models which describe the horizon fluctuations. In this 
paper, we will propose a different model, in which quantized linear perturbations 
of the gravitational field act as the source of the underlying metric fluctuations. 
Our analysis will be based on the formalism for the study of lightcone fluctuations 
proposed in Ref. ]4], and further developed in Ref. J5|. 

The necessary formalism will be reviewed and extended in Sect. ^. It will be 
applied to the case of cosmological horizons in Sect. [3], and to black hole horizons in 
Sect. [|. Our results will be summarized and discussed in Sect. [5| We will also give a 
critical assessment of the previous attempts 0, [| to estimate the horizon fluctuations. 

2 Basic Formalism 

In Ref. @], henceforth I, a model of lightcone fluctuations on a flat background was 
developed. It was assumed that the quantized gravitational field is in a squeezed 
vacuum state. This is the natural quantum state for gravitons produced by quantum 
particle creation processes, as for example in the early universe. Here we wish to 
generalize this formalism to the case of curved background spacetimes. Consider an 
arbitrary background metric gfij with a linear perturbation h^ u , so the spacetime 
metric is |J 

ds 2 = (g$ + h^dx^dx" . (1) 

For any pair of spacetime points x and x', let a(x,x') be one half of the squared 
geodesic separation in the full metric, and a (x,x') be the corresponding quantity in 
the background metric. We can expand cr(x, x') in powers of as 

a = a + 0i + cr 2 H , (2) 
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where o\ is first order in h^ v , ect. We now suppose that the linearized perturbation 
hfj, v is quantized, and that the quantum state is a "vacuum" state in the sense 
that we can decompose h^ v into positive and negative frequency parts 7t+ and h~, 
respectively, such that 

^> = 0, (Vl^ = 0. (3) 

It follows immediately that 

( V) = (4) 

in state In general, however, ((h^ u ) 2 ) ^ 0, where the expectation value is under- 
stood to be suitably renormalized. This reflects the quantum metric fluctuations. 

We now wish to average the retarded Green's function, G ret (x,x'), for a massless 
field over the metric fluctuations. In a curved spacetime, G ret (x,x') can be nonzero 
inside the future lightcone as a result of backscattering off of the spacetime curvature. 
However, its asymptotic form near the lightcone is the same as in flat spacetime: 

G ret (x,x')~ d Jtz±} 5 (cr), a^O. (5) 

47T 

We will ignore the backscattered portion, and average this delta-function term over 
the fluctuations, following the method of I. The result is 

(GM) = «p(-^) . (6) 

The effect of the averaging has been to replace the delta-function by a Gaussian with 
a finite width determined by the magnitude of the quantity (erf), which is the measure 
of the lightcone fluctuations. 

The operational meaning of the smeared lightcone can be understood by con- 
sidering a source and a detector of photons. If we ignore the finite sizes of photon 
wavepackets, then in the absence of lightcone fluctuations, all photons should tra- 
verse the interval between the source and the detector in the same amount of time. 
The effect of the lightcone fluctuations is to cause some photons to travel slower 
than the classical light speed, and others to travel faster. The Gaussian function in 
Eq. (|B|) is symmetrical about the classical lightcone, a = 0, so the quantum lightcone 
fluctuations are equally likely to produce a time advance as a time delay. 

In order to find the magnitude of the lightcone fluctuations in a particular sit- 
uation, it is necessary to calculate o"o for the metric in question, as well as (af) in 
the appropriate quantum state. This enables one to find At, the mean time delay 
or advance (measured in a suitable reference frame). This is an ensemble averaged 
quantity, not necessarily the expected variation in flight time of two photons emitted 
in rapid succession. To find the latter quantity, one must examine a correlation func- 
tion. This is the topic of Ref. 0. In the present paper, we will not be concerned with 
correlation functions, and will use Eq. to estimate the magnitude of the horizon 
fluctuations. 
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We may find a general expression for (erf), which is the curved space generalization 
of the result obtained in I. Let us first consider timelike geodesies. If we adopt a 
timelike metric, then in Eq. (|1|) we have that ds 2 > 0. Let w M = dx^/dr be the tangent 
to the geodesic and r be the proper time. We will define (af) by integrating along 
the unperturbed geodesic, in which case w M is normalized to unity in the background 
metric: 

$ } «V = 1 . (7) 
The geodesic interval in the unperturbed metric is given by 

^o = ^(Ar) 2 , (8) 

where At is the proper time elapsed along the geodesic. We have 

^ = yjl + V u»u v « 1 + hi^ u»u v , (9) 

and hence the geodesic length between a pair of points in the perturbed metric is 
As = At + Asi, where 

^ = \j dT Kv^u u . (io) 

Thus 

a = i(A S ) 2 = 1 -{At) 2 + ArAsi + 0{h 2 ) , (11) 
and hence o\ = AtAsi. If we average over the metric fluctuations, the result is 

= T^o J dT X dT 2 u^u\u p 2 u 2 (h^(x 1 )h pa (x 2 )) , (12) 

where Ui = dx^/dri and u 2 = dx^ jdT 2 . An analogous expression holds for the case of 
a spacelike geodesic, in which the integrations are over the proper length parameter 
of the geodesic: 

(erf) = --<To J d\ 1 d\ 2 utu±u 2 °u 2 T (h^x^hpvfa)) , (13) 

where now = dx^/dXi is the tangent to the geodesic, and A is the proper length. 
Here we have a = — |(AA) 2 . 

As noted previously, the quantity (a 2 ) is formally divergent and needs to be renor- 
malized. This may be done by defining the graviton two-point function, {h pu {xi)h pa {x 2 )) 
using, for example, the Hadamard renormalization scheme proposed by Brown and 
Ottewill 0]. These authors give a detailed prescription for expanding the singular, 
state-independent parts of the scalar and vector two-point functions in an arbitrary 
curved spacetime. Hadamard renormalization consists of subtracting this expansion 
from a given two-point function. This procedure seems not to have been developed 
in detail for the graviton two-point function, but there seems to be no barrier in prin- 
ciple to doing so. Allen et al M have applied the Hadamard renormalization method 
to the graviton two-point function in the vicinity of a cosmic string. In this paper, 
we will be content with simple approximations or order of magnitude estimates, and 
will not require the full renormalization machinery. 
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3 Cosmological Horizons 



Consider a spatially flat Robert son- Walker universe, for which the metric may be 
written as 

ds 2 = a 2 {r]){dr] 2 - rfx 2 ) . (14) 

In general, this spacetime has "particle horizons" associated with the comoving ob- 
servers across which other observers may appear or disappear. In the case of a 
radiation or matter dominated universe with an initial singularity (a oc rj or a oc i] 2 , 
respectively), a given observer at a given time has not yet received any light signals 
from distant observers, who are said to be outside of the first observer's horizon. In 
the case of de Sitter space (a oc r/^ 1 ), a given observer eventually ceases to receive 
signals from other comoving observers, and views them as having moved outside of 
the horizon. Clearly, these cosmological horizons are observer dependent in a way 
that black hole event horizons are not, and are basically the past lightcone of a given 
observer at a given time. Nonetheless, it will be of interest to estimate the magnitude 
of the quantum fluctuation of these horizons in various models. 

We must first study timelike and spacelike geodesies in the limits in which these 
approach null geodesies over some interval. Because the lightcone fluctuations are 
symmetrical, we may focus our attention on the timelike case. The geodesic equations 



for a timelike observer moving in the x-direction in the metric of Eq. (14) may be 
expressed as 

a2 dx = 1 
dr 



and 

d 2 ri a! f dr)\ a! . . 

where r is the proper time along the geodesic, a 1 = da/dt], and a is a constant. In 
the flat space limit (a = 1), we find from Eq. (|15|) that a = (1 — v 2 )/(2v 2 ) where v is 
the magnitude of the three-velocity. Thus in the null limit, a — > 0. The geodesic for 
a particle which starts at r\ = r/o at x = may be expressed as 

x = r] - rjo - f(v,Vo) , (17) 

where / — > in the null limit. For nearly null geodesies, we may assume |/| 1. 
To first order in /, the solution of Eqs. ( |i~5l) and (^6|), corresponding to a properly 
normalized four-velocity, is 

f(v,Vo) = a [ V a 2 ( V ')d V ' . (18) 



From Eq. (|i4[), we can write 



ds = \jl— (— ^ a dr] rs \p2a a 2 dr/ , (19) 
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from which we find 

At = V2a / a 2 (r))dr] , (20) 

and 

^o = ^(Ar) 2 . (21) 

Gravitons propagating on a Robert son- Walker background may be quantized in 
the transverse, trace- free gauge, which eliminates all of the gauge freedom ||. Only 
the purely spatial components of are nonzero, and they each satisfy the wave 
equation for a massless, minimally coupled scalar field in this metric. Thus the 
graviton two-point function can be expressed in terms of the scalar two-point function, 
(ip(x)ip(x')}, as 

1 3 3 

(hij( x )h k i(x')} = --a 2 (r?)a 2 (?/)(%5 M - -5 ik 5ji - -SuSjk) (<p(x)ip(x')} . (22) 

We may use this result to write 

= [ m d V [* d V ' (<p(xMx')) , (23) 

Note that in this expression, we can replace the unsymmetrized two point function, 
{(p(x)ip(x r )), by the symmetrized form (the Hadamard function) 

G(x,x') = ±{<p{x)<p(J) • ( 24 ) 

The latter function is real, and is hence more convenient. To proceed further, we 
must designate the quantum state of the gravitons. In the following two subsections, 
some particular examples will be considered. 



3.1 Gravitons in a Radiation-Dominated Universe 

A radiation-dominated universe, for which 

a(r])=a r], (25) 

is presumably a reasonably good description of a significant fraction of the history 
of our universe. Let us consider a thermal bath of gravitons in such a universe, for 
which the temperature is always high compared to the scale set by the local radius 
of curvature, i.e., the thermal wavelength is much less than the horizon size. In this 
case, the minimally coupled scalar field two point function is approximately equal to 
that for the conformally coupled field, G cc (x,x'). However, the latter is conformally 
related to the flat space Hadamard function, Gq(x, x'), so we have 

G(x, x) ~ G cc (x, x') = a~ 1 (r^) a^ 1 ^') G (x, x') . (26) 
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We now need the flat space renormalized thermal Green's function on the lightcone. 
This was calculated in Appendix A of Ref. [El], where it was shown that in the high 
temperature limit, this function is given by 



G (x,x')~— - /?«p, (27) 

where (3 is the inverse temperature and p — |x — x'|. We may now use Eqs. (|23|) - 
( p7|) to write 

n On H 1 . 1 

of = / dri / d7 7 / , r . (28) 

Unfortunately, this integral diverges because of the singularity of the integrand at 
7] = rj . This is due to the fact that Eq. (|27j) is not valid for small p. We can remedy 
this by excluding the range \rj — rj'\ < e, where e is a cutoff which will be taken to be 
of order f3. Thus, the relevant integral is, in the limit of rji ^> r] ; 

, 9v on r?i , / r r ni \ 1 a o ln (^o/e) + i / s 

(of) = °— -z \ dr) + )dr}' — , r w — ^ ^ . (29) 

487r a (3 % Jrto \Jr] Je J VV'lV ~ V \ 247r a (3 7i clq 

In this same limit, one finds 

ao^-aalri^. (30) 

We now wish to define a c as that value of a for which the argument of the 
exponential in Eq. (Q) is unity, that is, 

^ 1 . (31) 



Thus a c describes a geodesic whose deviation from the lightcone characterizes the 
fluctuations. From Eqs. (^), (0), and (|31|), we find that 



]n(r j0 /e) + l j 

a c = Wn \ ^ ■ ( 32 ) 

We need to define a physical measure of the magnitude of the lightcone fluctuations. 
This may be taken to be the mean time delay or advance, Arj, for a photon emitted 
at 1] Q and detected at rjx. Equivalently, we can think of Arj as the characteristic 
interval around r/o within which photons could be emitted and still reach a detector 
at a coordinate distance of Ax = rji — r] at time rji. (See Figure 1.) From Eqs. (|l7|) 
and fll8|) , we have that Arj is related to a c by 

rvi 1 

A?7 = f(vi,Vo) = oc c / a 2 (r]')dr]' ^ -a c a 2 r]f . (33) 
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Ax 




x 



Source Detector 



Figure 1: 

A photon is received by a detector at conformal time rji. In the absence of metric 
fluctuations, it has traveled along the classical lightcone (dashed line) from a source 
at a coordinate distance |Ax| = rji — rjo, and was emitted at conformal time 770- In 
the presence of metric fluctuations, it could have been emitted a characteristic time 
A77 before or after 770, and traveled along a mean trajectory which is either timelike 
or spacelike, respectively (dotted lines). 

It is perhaps more convenient to express this time delay or advance as a coordinate 
time interval, At = 0(770) At?, which is given by 



/2t \i ^[HVoM + l 

At 



do / 6V7T/3 



(34) 



where to = |°o^o * s the coordinate time at r] . 

We may interpret this formula by noting that in an expanding universe, f3~ l is a 
coordinate temperature, not in general the physical temperature. It is, however, the 
physical temperature at a time at which a — 1. Let us take that time to be t , the 
time of emission, at which time the physical temperature is T . This leads to 

At V6 r— - ; ; /T \ft p 



^/NWOTI (£)(£)• (35) 



t V \J. p 



where T p is the Planck temperature, and t p is the Planck time. The logarithmic 
factor can be taken to be of order one, so we see that if T = T p and to —t p , then we 
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have At /t ~ 1. Otherwise, with T <C T p and t > t p , we have At/t <C 1, and the 
fractional lightcone fluctuations are small. This is what we should perhaps expect; 
a bath of gravitons with the Planck temperature at the Planck time (which would 
correspond to a few degrees Kelvin today) results in large horizon fluctuations, but 
otherwise the fluctuations are small at much lower temperatures. 



3.2 Gravitons in de Sitter Space 

If we represent de Sitter space as a spatially flat Robertson- Walker metric, then the 
metric is of the form of Eq. (]TJ|) with a(r]) = —H/rj = H/\rj\, where if is a constant, 
and — oo < r\ < 0. These coordinates cover one-half of the full de Sitter spacetime, 
but that is sufficient for our purposes. Gravitons are again represented as a pair of 
massless, minimally coupled scalar fields. However, in this case there is a subtlety in 
that there is no de Sitter invariant vacuum state which is free of infrared divergences 



10 , 11 , |12 |. The Hadamard function may be represented as [12 



G(x,x') = -^Re J d 3 kMv)r k (v')e i( -^ x , (36) 

where the time part of the mode function is expressible in terms of Hankel functions 

as 

Mv) = ^~H Ml \ Cl HP{k V ) + c 2 M 2) (fcr/)l . (37) 

2 1 2 2 J 

Here c\ and c 2 are functions of k which are required to have the following properties: 

ci — > 0, c 2 — > 1, as k — > oo , (38) 

and 

|ci + c 2 |->0, as k^O. (39) 



These functions define the quantum state in question. Equation ( |38"D insures that 
the very high frequency modes are free of particles, whereas Eq. (^) is the condition 
that the state be free of infrared divergences. 

The Hadamard function given by Eqs. (^) and ([IT]) is the unrenormalized function 
which is singular on the lightcone. In principle, one should evaluate the integral for 
the given choice of C\ and c 2 , and then extract the state- independent singular terms 
to obtain the renormalized Hadamard function. Unfortunately, this would be very 
difficult to do explicitly. However, in the late time limits (r/ —>■ or rj' — > 0) it is 
possible to give some approximate forms for the renormalized function. It was shown 



by several authors that the coincidence limit grows logarithmically |IT| [12], 0, [14 



G(x,x)~-^\nH\ri\, (40) 

as 7] = rj' — > 0. The fact that this asymptotic form is state-independent may be 
understood as a consequence of the exponential expansion having redshifted away 
any memory of the quantum state. 
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We may use a similar procedure to investigate the behavior as 77 — > with 7/ 
fixed. If we insert Eq. fl3~?|) into Eq. ([36]) and then change the variable of integration 

\v 



to q 

G(x, x') 



rj'\ k, we have 
H 2 \t]t]'\% 



327T 2 1 77 — T}' 1 3 



Re / d 3 q 





rf 






V 




rf\ 



2 



1*7 1 



1 77 — 7]' 





V 






V 




v'\ 





rf 






V 




T)'\ 



Q 



3 ivq 



(41) 



where v = (x— x')|r7— Now assume that 77 is sufficiently small that the dominant 
contribution to the integral comes from values of q for which the magnitude of the 
arguments of the Hankel functions in the first factor are small compared to unity In 
this case, we can use the small argument forms 



'2 _3 
X 2 
7T 



< x < 1, 



in this factor and find 
V2H 2 



G(x, x') 



32tt 5 /2 



Re / d 3 qg"i {c x + c 2 )[c\Hf (-q) + c* 2 Hf(-q)} e 



r(l)/ 

2 



,'v-q 



(42) 



(43) 



Note that the dependence upon both 77 and 77' has dropped out. The integral in 
Eq. (|43D has a logarithmic ultraviolet divergence on the lightcone, but the leading 
quadratic divergence has disappeared. A similar disappearance of divergent parts oc- 
curs in the derivation of Eq. ([40]). (See, for example, Ref. fl~2|| .) The renormalization 
of this logarithmic divergence requires a subtraction of a term proportional to the 
scalar curvature R = 12H 2 . We expect the result to be a constant which is of order 
of H 2 . By symmetry, we obtain the same result if 77' — ^ with 77 fixed, so that 



G(x, x) 



77 — > or 77' — > 



(44) 



The integral in Eq. Q2"3] ) requires us to know G(x, x') in the square illustrated in Fig. 
2, which also illustrates the regions in which either Eq. ( ^0|) or Eq. ([44]) is applicable. 
For the purposes of obtaining an order of magnitude estimate for (a 2 ), we will assume 
that G(x, x') is of the form 



G(x,x') ss H 2 F(rj,ri') 



(45) 



throughout this region, where F is either a constant of order unity or else a logarithmic 
function which will contribute multiplicative constants of order unity to the integral 
in Eq. (p3|). 

The result of the evaluation of this integral is 



(*?> 



H 2 a 

K (t?i 

a 



(46) 
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^1 
-1 




1i 



Figure 2: 

The domain of integration for (af) is the interior of the square. The integrand, 
G(x,x f ) is known in the shaded regions to be approximately H 2 , and in the 
crosshatched region to be approximately — (H 2 \nH\r]\)/4:7Y 2 . 
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where k is a constant of order unity. From Eq (p0|) we find 

g(7yi - T] ) 2 

ao = ^¥F- (47) 

The value of a which is characteristic of the lightcone fluctuations is 

a c « V2kH 3 T]iT] . (48) 

If we follow the line of reasoning in the previous subsection, we find that the charac- 
teristic time associated with the de Sitter horizon fluctuations is 

At^V2n. (49) 

This is of order unity, and hence again the horizon fluctuations are of Planck dimen- 
sions. 

We may use At to find the frequency fluctuations observed at rjx from a constant 
frequency source. Suppose that the source is emitting photons with a constant fre- 
quency uq. In the absence of metric fluctuations, the photons will be detected at 
frequency v = vq a(r]o) /cl(t]i)- The effect of the metric fluctuations is equivalent to 
a drift in the source frequency whose magnitude is Au = v$ At. Consequently, the 
fractional variation of frequency at the detector is 

Av 

— = u At. (50) 

However, this is an ensemble averaged frequency variation, not necessarily the drift 
in frequency that would be observed in any one trial. The reason for this is that 
pulses emitted close together in time tend to have correlated time delays or advances 
f|. Thus, Eq. (|50| ) should be interpreted as giving an upper bound in the frequency 
drift seen by the detector. 



4 Black Hole Horizons 

Here we wish to discuss the fluctuations of the event horizon of a Schwarzschild black 
hole, for which the metric is 

ds * = (i_^ dt 2_ (i_ MlyV-r 2 (^ 2 + sin 2 ^# 2 ). (51) 

Timelike radial geodesies in this metric satisfy [[U| 

(!)' = #-C(r) (52, 

and 

?- = E/C, (53) 
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where C(r) = l — 2M/r and E is a constant of the motion which is equal to the energy 
per unit rest mass of the particle, as measured at infinity. From these relations, one 
may show that 

and that 



where 



dtj V E 

dr* ~ ' ' ' 

From Eq. (p3p, we see that the proper time elapsed along a segment of a geodesic is 

At = E- 1 I 1 C(r)dt, (57) 

Jto 

where r is understood to be a function of t along the geodesic. We are primarily 
interested in the case of nearly null outgoing geodesies, for which E — E ^> 1. In 
this limit 

(IT 

- « C(r) . (58) 

Thus, |Ar| w |Ar|/E and 

where Ar is the radial coordinate interval traversed. An analogous treatment may be 
given for spacelike geodesies. The constant of the motion E no longer has a simple 
physical interpretation, but we can express a in this 

(Ar) 2 , s 

2E 2 V > 

We now turn to the task of estimating (af) near the horizon of a Schwarzschild 
black hole. From Eq. ( [T2]) and the fact that 

-r^E Q (6i) 

GST 

for nearly null outgoing timelike geodesies, we have that 

~ ^°" ^ 2 / rfri rfr2 U \ U \ U 2 U 2 (h^u{xi)hp a (x 2 )) ■ (62) 

If we are interested in a black hole radiating into empty space, then the relevant 
quantum state for the quantized graviton field is the Unruh state. It would be a 
rather formidable task to explicitly compute the renormalized graviton two point 
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function in the state. Instead, we will content ourselves with an order of magnitude 
estimate. First we must choose a convenient gauge. Again we wish to impose the 
transverse, tracefree gauge, which eliminates all gauge freedom. Because all of the 
modes of the graviton field are propagating waves which either originate at Z~ or 
reach X + , we can impose the requirement that these modes satisfy the flat space 
transverse, tracefree gauge condition at r = oo. 

We now make the assumption that in this gauge, the renormalized two point 
function measured in the frame of an infalling observer who starts from infinity at 
rest can be estimated by dimensional considerations. Near the horizon at r = 2M, 
the geometry and the quantum state are characterized by a single scale, M. If we 
were to reinstate explicit factors of Newton's constant, G, then h^ u oc oc m p 
where m p is the Planck mass. However, is dimensionless in any set of units, 
so our assumption tells us that the two point function should be proportional to 
rrip/M 2 . However, the actual values of the components of this bitensor depend upon 
the choice of frame. Our assumption is that infalling observers with E = E^ ~ 1 
should be regarded as preferred in the sense that they do not introduce any very 
large or very small dimensionless redshift or blueshift factors. Let be the four 
velocity of such an observer. Our assumption may be expressed as 



2 

vat 



(M*i)M*a)> « j£ ( 63 ) 
near the horizon. The components of the infalling observer's four velocity are 



v 



C-\ v r = -VT^C « -1 , (64) 



and those for an outgoing observer are 

u l = E C'\ u r = y/E*-C&E . (65) 
Thus | | = E Q | v ^ | and we can write 

2 

~ Tn 

«f«M«? (M^OM 2 *)) ~ E oj^- ( 66 ) 



Our basic assumption receives some support from the work of York [[15] who esti- 
mates the magnitude of the quantum fluctuations of the lowest modes of vibration of 
a Schwarzschild black hole. He treats these modes as quantum mechanical harmonic 
oscillators, and calculates their root-mean-square fluctuation amplitudes. The am- 
plitudes of the first few modes yields a result consistent with Eqs. ( |63D or (|66|). Of 
course, this is heuristic support, and by no means a proof of our assumption. A full 
proof would require one to sum over an infinite number of degrees of freedom, and 
then extract any ultraviolet divergent parts. 

The graviton two point function in this approximation is a constant in the vicinity 
of the horizon. It must also fall off to zero at large distances from the black hole. Thus 
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the integral in Eq. (BS1) gets its dominant contribution over an interval in r of the 
order of M, regardless of the upper limit of the integration. In any case, we can stop 
the integration at a maximum value of r which is just a few times M. Whether the 
outgoing photons emitted in the vicinity of the horizon are detected at r = 4M or at 
a much larger value of r has little effect on the discussion of the horizon fluctuations. 
Thus we may let 

W)««r„B^. (67) 

In analogy to the discussion in the previous section, we wish to define a characteristic 
value E c , which is the value of E Q at which Eq. (^l|) holds. From Eqs. (|59|) and (|67l) 
we find 



E c » VM . (68) 

We may find the associated time delay or advance, At, from Eq. ([55]), which tells us 
that when E ^> 1 

C 1 
dt « dr* + - dr* = dr* H =- dr . (69) 

A radial null geodesic in the classical background geometry covers an r* distance of 
Ar* in a coordinate time At = Ar*. The second term on the right hand side of 
the above equation tells us the extra amount of time required by a timelike particle. 
Analogous expressions hold for spacelike geodesies, and yield the same magnitude of 
time variation. Thus we are led to an expression for the characteristic time delay or 
advance due to horizon fluctuations: 

Ar 

As discussed above, we can take Ar to be of order M, although we might also want 
to consider the possibility of taking it to be much smaller. Thus let 

Ar = 7 M , (71) 

where 7 is a constant of the order of or less than unity. Now we have 

At « 7 , (72) 

so the time delay, measured in coordinate time, is of Planck dimensions. However, 
a more physical measure is obtained by expressing this time interval in terms of the 
proper time of a local observer. Let the photons be emitted at r = r = 2M(1 + e), 
with e C 1, and let Co = C(r'o) ~ e. The time interval in the frame of a static 
(nongeodesic) observer at rest at r = r is 

AT.nyy/c , (73) 
that in the frame of an infalling observer with E — 1 is 

Ar i ^ 7 C ^7e, (74) 
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and that in the frame of an outgoing geodesic observer with E = 1 is 

Ar « 7 . (75) 

One might regard Arj, the characteristic time as measured by an infalling observer, 
to be the best measure of the magnitude of the horizon fluctuations. Such an observer 
can cross and continue beyond the classical event horizon at r = 2M. Suppose that 
an outgoing photon emitted by this observer reaches infinity. An observer at infinity 
who detects this photon and who is unaware of the lightcone fluctuations might trace 
the history of this photon backwards in the classical Schwarzschild geometry and infer 
that it was emitted at a proper time of tq on the infalling observer's worldline. In fact, 
it could have been emitted anywhere in a band of width Ar^ centered around To- (See 
Fig. 3.) The remarkable feature of the result Eq. ([71p is that Ar,, — > as r — » th, 
the proper time at which the infalling observer reaches r = 2M. In the cosmological 
models discussed in Sect. |3|, the fluctuation in emission time was typically of the order 
of the Planck time. In the black hole case, the horizon fluctuations are more strongly 
suppressed. Note that the proper time required for the infalling observer to pass from 
r = r = 2M(1 + e) to r = 2M is T ~ 2eM. This is always large compared to Arj 
for large black holes: 

Arj m„ , , 

— - « 7^ . 76 
T M K J 

Thus the only outgoing photons which manage to cross the classical horizon are part 

of an extreme tail of a Gaussian distribution. 

As in Sect. |37|, we may express the time delay or advance in terms of the variation 

in frequency seen by the observer at infinity. In the black hole case, the analog of 

Eq. (!§ is 

— = uoAn w z/ 7e- (77) 
v 

Thus as the source approaches r = 2M , the fractional variation in frequency observed 
at infinity goes to zero, and the observed frequency approaches that predicted by 
classical relativity. 

Let us now turn to the question of whether horizon fluctuations are capable of 
invalidating the semiclassical derivation of the Hawking effect. First let us recall 
the essential features of this derivation, as given in Hawking's original paper ||17|| . 
Consider the spacetime of a black hole formed by gravitational collapse (Figure 4). 
The null ray which forms the future horizon leaves Z~ at advanced time v — v q. The 
modes into which the outgoing thermal radiation will be created leave Z~ at values 
of v slightly less than t> , pass through the collapsing body, and reach 1 + as outgoing 
rays, on which the retarded time u is constant. Hawking shows that the relation 
between the values of v and of u is 

u=-4Mln(^^), (78) 

where A is a constant. Thus u — ► oo as v — ► vq. As seen by an observer at infinity, 
these outgoing rays must hover extremely close to the horizon for a very long time. 
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Inf ailing 
Observer 



Figure 3: 

An observer falling across the future horizon, H + , of a black hole emits photons which 
reach J + . In the presence of metric fluctuations, these photons need not follow the 
classical lightcone (solid line), but rather may follow timelike or spacelike paths in the 
background geometry (dotted lines). The characteristic variation in emission time, 
as measured in the frame of the infalling observer, of photons which reach X + at the 
same point is Ar^. 
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Figure 4: 



The spacetime for a black hole formed by gravitational collapse. The shaded region 
is the interior of the collapsing star. A null ray which leaves X~ with advanced time 
Vo becomes the future horizon, H + . A ray which leaves at an earlier time v passes 
through the collapsing body and reaches X + at retarded time u. The dashed line is 
the worldline of a observer who falls into the black hole after its formation. 

If one starts with a black hole with a mass M large compared to the Planck mass, 
the semiclassical description should hold for the time required for the black hole to 
lose most of its original mass. Let 



be this characteristic evaporation time. The basic problem posed by the horizon 
fluctuations is that they may cause an outgoing ray to either fall back into the black 
hole, or else to prematurely escape. In either case, the semiclassical picture of black 
hole radiance would need to be modified at times less than t evap . 

At a large distance from the black hole, u = t — r* m t — r. If the observer at 
"infinity" is at a fixed value of r (e.g. 100M), then u ~ t for most of the black hole's 
lifetime. Thus, in order not to invalidate the semiclassical treatment, outgoing rays 




(79) 
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with u < u max = t evap need to be uninfluenced by the horizon fluctuations. In order 
to investigate this question, let us consider an infalling observer with E = E\ = 1. 
From Eq. fl54|), we have that near r = 2M 

I " ~ c " < 80 > 

and hence 

^ = #-^«-2CT 1 . (81) 
dr dr dr 

This equation may be integrated to yield 

u(r) = -4Mln( — —J, (82) 

where A' is a constant. This relation tells us the value of r at which the infalling 
observer crosses a given constant u line. The constant A' is determined by which 
infalling observer we consider. Here we are interested in observers who fall into the 
black hole not long after its formation, and we can set A' rs M. Let r c be the value 
of r at which this observer crosses the u = u max line, given by 

r c — 2M = M e ~ Umax/4M . (83) 

Near the horizon, Eq. fl52|) tells us that dr/dr w —1 along the worldline of the infalling 
observer. Thus the proper time required for this observer to cross from u = u max to 
the classical horizon at r = 2M is 

5t w r c - 2M = M e ~ u ™-/ 4M « M e- M2/m ? . (84) 

We should compare this quantity with Ar^, where C is evaluated at r = r c , so 
C = I e -«— / 4M . Thus, 

Ar ? = ^5r, (85) 

2M v ; 

and hence so long asM> m p , Arj <ti St. From this result, we conclude that the 
horizon fluctuations do not invalidate the semiclassical derivation of the Hawking 
effect until the black hole's mass approaches the Planck mass. This is the point at 
which we would expect the semiclassical treatment to fail. 

The presence of frequencies far above the Planck scale, in the form of the modes 
leaving X~, has concerned numerous authors. There have been suggestions that one 
might be able derive the Hawking effect in a way that transplanckian frequencies do 
not arise, using some form of "mode regeneration" [|T8j, So far, it has not been 
possible to implement these suggestions in detail. As seen from the our analysis of 
horizon fluctuations, the semiclassical treatment is remarkably robust. 
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5 Summary and Conclusions 



In the preceeding sections, we have analyzed the horizon fluctuation problem using a 
formalism which takes account of the effects of quantized linear perturbations of the 
gravitational field upon lightcones. In the case of the cosmological models treated in 
Sect. ^ the resulting horizon fluctuations were found to be of Planck dimensions for 
both de Sitter space and a radiation filled universe with a Planck density of gravitons 
at the Planck time. These fluctuations are measured as fluctuations in the time of 
emission of a photon as measured in the frame of a comoving observer. The order of 
magnitude of the results is what one might have guessed before doing the calculation. 

In the case of black hole horizon fluctuations, the results are somewhat more 
subtle. Whether the time scale which characterizes the horizon fluctuations (the 
time delay or advance) is of Planck dimensions or not depends crucially upon the 
frame of reference. It is indeed of Planck dimensions as measured by an observer 
at infinity. However, as measured by an infalling observer, this time is much less 
than the Planck scale, and vanishes as the infalling observer approaches the classical 
event horizon at r = 2M. We further found that this suppression of the horizon 
fluctuations is exactly what is needed to preserve Hawking's semiclassical derivation 
of black hole radiance for black holes of mass large compared to the Planck mass. 

Our result seems to conflict with the arguments of Sorkin |2|] and of Casher et al 
H . These authors claim that the horizon fluctuations are much larger than found in 
the present manuscript. It should be noted, however, that the physical mechanisms 
being postulated in Refs. [fj] and are quite different from that of the present 
paper. Furthermore, in our opinion, the physical basis of both of these calculations 
seems to be open to question. Casher et al obtain large gravitational perturbations 
of the horizon by postulating an "atmosphere" of particles near the horizon in large 
angular momentum modes. This arises by decomposing the physical quantum state 
of an evaporating black hole (the Unruh vacuum) into two pieces which separately 
have divergent stress tensors on the horizon, the contribution from the Boulware vac- 
uum state and a term which these authors call the "atmosphere" of particles. The 
large stress tensor fluctuations arise in the Casher et al analysis when this "atmo- 
sphere" undergoes thermal fluctuations. Our objection to this procedure is that the 
fluctuations of the Boulware vacuum energy density are not being considered. The 
splitting of the finite Unruh vacuum energy density into two singular parts seems 
rather artificial. If one chooses such a splitting, then care must be taken to prove 
that fluctuations in one part are not cancelled by correlated fluctuations in the other 
part. Casher et al have not done this. 

Sorkin 0] uses a Newtonian treatment to estimate the gravitational field of a mass 
fluctuation near the horizon, and its effects on the Schwarzschild geometry. One 
can certainly question whether a Newtonian analysis can be trusted in black hole 
physics. However, our primary objection to Sorkin's treatment is that the dominant 
contribution to the horizon fluctuations comes from modes whose wavelength is very 
small compared to the size of the black hole. The same line of reasoning would seem 



20 



to lead to large stress tensor fluctuations, and hence large lightcone fluctuations, 
in all spacetimes including flat spacetime. In our view, a more reasonable result is 
one in which significant fluctuations arise only on scales characterized either by the 
spacetime geometry, or else the chosen quantum state. An approach to defining stress 
tensor fluctuations on a flat background which has this property was given in Ref. 
|20|| . Here the stress tensor fluctuations are defined in terms of products of operators 
which are normal ordered with respect to the Minkowski vacuum state. 

Recently, the fluctuations of the Hawking flux, as measured in the asymptotic 
region, have been computed |2T| by a similar approach. It was found that this flux 
undergoes fluctuations of the same order as its average value over time scales of the 
order of M. This average flux is of order M~ 2 , so the characteristic associated black 
hole mass fluctuation is of order M~ l . The corresponding metric fluctuation near the 
horizon is then of order 5h « M~ 2 . For macroscopic black holes, this is much smaller 
than the metric fluctuations due to the quantized linear perturbation, estimated in 
Eq. ( |63"D to be of order M~ l . This analysis does not rule out the possibility of much 
larger stress tensor fluctuations in the vacuum energy near the horizon. However, the 
diagonal and off-diagonal components of the expectation value of the stress tensor in 



the Unruh state near the horizon are of the same order |[22|| . It is thus plausible that 
the fluctuations in these various components near the horizon are also of the same 
order. If so, then the effects of quantized linear perturbations of the gravitational 
field dominate over those of stress tensor fluctuations. 

It must be emphasized that all of the conclusions obtained in the present manuscript 
are in the context of a model of linearized quantum gravity. Furthermore, much of 
our discussion is of a heuristic, order-of-magnitude nature. If the basic picture of 
horizon fluctuations which we have drawn is correct, much work remains to be done 
to make the picture more precise. 
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